The Plant Journal

Whether Photosystem I (PSI) core subunits accumulate prior to light exposure in developing angiosperm seedlings remains unresolved, with conflicting reports across species. Here, we investigated the presence and membrane colocalization of the PSI core subunit PsaA in etioplasts of dark-grown angiosperms representing dicot and monocot species. Immunoblotting showed that PsaA accumulates in etioplasts of all three dicot species examined (pea, Arabidopsis, and runner bean), whereas in the monocot oat it was detected only after prolonged etiolation, at substantially lower levels and with an anomalously high apparent molecular weight. Blue-native PAGE analysis reveals that a fraction of PsaA co-migrates with LPOR, PsaB, FNR, and chlorophyll synthase, suggesting co-localization within a shared membrane microdomain rather than stable complex formation. The thylakoid insertase Alb3 was more abundant in dicot etioplasts, consistent with a potential role in the early integration of PsaA into the membrane. Upon illumination, pea reached PSI functionality faster than oat, with P700 oxidation detectable 30 min earlier, linking the dark accumulation of PsaA to an accelerated photosynthetic onset. These findings demonstrate light-independent accumulation of a PSI core subunit in a species-dependent manner and point to early steps in PSI biogenesis that precede full photosynthetic complex assembly. Highlight Contrary to prevailing models, a Photosystem I core subunit PsaA accumulates in dark-grown angiosperm seedlings before light exposure, revealing light-independent early steps in photosynthetic complex biogenesis.

The TETRATRICOPEPTIDE THIOREDOXIN-LIKE (TTL) proteins are a plant-specific family proposed to function as peripheral membrane proteins that contribute to abiotic stress tolerance in Arabidopsis, likely by maintaining cell wall integrity through brassinosteroid signaling. Previously, we identified a TTL gene that was differentially regulated at the translational level in nitrogen-fixing soybean plants under water deficit (WD) conditions. This finding prompted the characterization of the soybean TTL gene family. Using the Glycine max v4.0 proteome, we identified ten TTL homologs (GmTTL1-GmTTL10), which are unevenly distributed across five chromosomes. Phylogenetic and structural analyses grouped these genes into three clades and revealed a highly conserved exon-intron organization. Likewise, GmTTL proteins display a conserved number and arrangement of TPR and TRXL motifs. To gain insights into their potential biological functions, we integrated co-expression and differential expression analyses. This approach identified a co-expression module enriched for translationally downregulated genes related to the Gene Ontology terms "cellular anatomical entity", "membrane", "cell periphery", "cell wall modification", "nitrate assimilation", and "cell wall organization or biogenesis". Protein-protein interaction network analysis of this specific subset of genes uncovered a novel GmTTL connection with two nitrate reductase enzymes in nitrogen-fixing plants subjected to WD, potentially linking the TTL gene family to new functions or roles. This study provides a framework for future functional studies of GmTTL proteins and their contribution to abiotic stress adaptation in soybean. Key MessageThis work presents the first functional characterization of TTLs proteins in legume species and highlights key processes that may link the TTL gene family to new functions or roles.

O_LIThe control of amount of nuclear proteins is fundamental in regulating plant gene expression, but the mechanisms of quantitative dynamics of the nuclear proteome are largely unstudied during adaptive responses to pathogens. C_LIO_LIHighly specific labeling, enrichment and measurement of the nuclear proteome was performed using TurboID LC-MS of Arabidopsis thaliana leaves treated with the pathogen-associated molecular pattern (PAMP), flg22, and/or cycloheximide. The chosen experimental approach allowed discrimination of the effects of translation, nuclear protein import, trafficking of preexisting proteins, derepression, and nuclear protein turn-over upon elicitation of basal immunity. C_LIO_LIThe highly specific, deep coverage of proteins in the nucleus makes this study a resource for anyone interested in plant nuclear proteome dynamics and defense. C_LIO_LIAround 2,000 nuclear proteins were repeatedly quantified, including more than 300 transcription factors or other proteins related to transcription. Several proteins with documented activity in endosomes were newly synthesized and imported into nuclei upon PAMP challenge, suggesting alternative nuclear functions in PAMP-triggered immunity (PTI). Circadian clock components, including the transcription factor, CIRCADIAN CLOCK ASSOCIATED 1 (CCA1)-HIKING EXPEDITION (CHE), were depleted upon PAMP challenge, suggesting a safeguard against untimely induction of systemic acquired resistance (SAR). C_LIO_LIBased on proteomic patterns, proteins moonlighting in the nucleus as well as trafficking and turn-over regulation of the proteome are common elements during plant immunity. C_LI

Fluorescence Lifetime Imaging Microscopy (FLIM) is becoming a key technique for live-cell multiplexing and label-free detection of endogenous fluorescence in animal systems. Its potential in plant biology, however remains largely unexploited, despite its integration into a number of commercial microscopy setups. Here, we build a systematic, subcellular FLIM reference library for a panel of genetically-encoded fluorophores. Lifetime imaging of different fluorescent reporters targeted to distinct organelles (nucleus, plasma membrane, endoplasmic reticulum, etc.) and subsequent analysis of the decay curves using different modes allowed us to simultaneously discriminate up to four spectrally overlapping fluorophores solely by lifetime differences in specific subcellular compartments. Remarkably, fluorophores with lifetimes differing by as little as 0.1 ns can be reliably discriminated using one of these modes, namely Phasor-based analysis. Moreover, we show that the same fluorophores exhibit compartment-specific lifetime shifts, enabling Phasor separation of identical tags residing in different organelles. Finally, we extended the Phasor approach to label-free imaging of endogenous plant fluorescence. Together, these results establish FLIM-Phasor as a versatile, multiplex-capable tool for plant cell biology, opening new avenues for imaging strategies that yield higher content information at both cellular and tissue-level resolution.

Live imaging of plant subcellular structures is key to deciphering the spatiotemporal bases of cellular processes, and their functional impact on growth and morphogenesis at various biological scales. Live imaging of plant cells essentially relies on expression of fluorescent markers labeling cells or subcellular structures of interest. Simultaneous multi-channel imaging of several markers is still not routine practice in plant cell biology, owing to issues linked to genetic or spectral compatibility of markers, differences in expression levels, silencing, toxicity, etc. Here we designed a three-color marker in Arabidopsis thaliana and Capsella rubella, enabling high-resolution live imaging of plant morphogenesis, including labeling of the cell membrane, the nucleus and the microtubule cytoskeleton. Detection of MT arrays involved the development of a MAP4-MBD-based microtubule marker optimized for plant cells. The three-color marker allows visualization of the three-dimensional organization and dynamics of plant microtubules within the intracellular space with unprecedented precision, in various organs including the root and shoot meristems, the leaf, anther, and gynoecium. Our results demonstrate the potential of such single-construct strategy for cell biology studies in plants.

Genes encoding intracellular nucleotide-binding site leucine-rich repeat (NBS-LRR) receptors represent the largest class of resistance (R) genes in plants, yet their evolutionary trajectory in trees remains poorly understood. Using high-quality long-read genome assemblies from eight Myrtaceae species, we identified 15,792 NBS-encoding genes with threefold variation in gene content across species. To investigate potential decoy domains for effector proteins of pathogens, we parsed the annotated NBS-encoding genes and determined 181 unique non-canonical domains. Notably, we observed frequent integration of Jacalin domains into TIR-NBS proteins. This novel gene family, named TNJ, is hypothesized to represent a new R gene class. Further exploration of TNJ sequence analysis shows up to seven Jacalin domains per protein and forming a monophyletic clade; however, AF3 modelling confirmed only six domains. Conserved residues in the TIR domain and functional motifs within the NB-ARC domain supports TNJs potential role in immune signalling. Most hypervariable sites and positively selected sites detected to be surface-exposed were clustered in the Jacalin region of TNJ, suggesting that surfaces of Jacalin domain may harbour residues determining pathogen recognition specificity. These findings support TNJ as a potentially new class of R genes in Myrtaceae with Jacalin as a replacement for LRR.

Plant immunity is, amongst others, mediated by receptor like proteins (RLPs), which are localized on the plasma membrane. RLPs recognize extracellular immunogenic patterns (ExIPs) originating from pathogens or derived from the host itself, which leads to extracellularly triggered immunity (ExTI). Cf proteins, which are well-known RLPs of tomato (Solanum lycopersicum) confer resistance against the fungal pathogen Fulvia fulva. Cf-9, Cf-4, Cf-2 and Cf-5 are well-known examples of Cf proteins, mediating recognition of the matching ExIPs Avr9, Avr4, Avr2 and Avr5, respectively, which are secreted effectors of F. fulva and trigger hypersensitive response (HR)-related cell death in tomato plants carrying these Cf proteins. Although all these Cf proteins confer proper resistance to the fungus, Cf-9 and Cf-4 trigger a stronger and faster cell death than Cf-5 and Cf-2. It is unknown whether these phenotypical differences arise from variations in the molecular mechanism of the cellular immune response that is initiated by the Cf proteins, and whether this phenotypic difference correlates with varying degrees in the intensity and timing of the triggered immune responses and robustness of the resistance. To try to answer these questions, in this study the immune responses triggered by Cf-4 and Cf-5 were compared. Cf-4 and Cf-5 share the same core upstream signaling components to trigger HR-related cell death in Nicotiana benthamiana. In tomato, both receptors induce rapid MAPK activation, which is more sustained for the Cf-5/Avr5 combination. Both Avr4 and Avr5 induce an apoplastic burst of reactive oxygen species (ROS), independently of the presence of their matching receptors, while remaining dependent on RBOHB for this ROS burst. Full transcriptome analysis at 3 and 7 hours after immune activation revealed a large shared set of differentially expressed genes, alongside qualitative and quantitative differences, with the Cf-5/Avr5 combination inducing a broader transcriptional reprogramming. Despite these differences, Cf-4 and Cf-5 confer a comparable level of resistance to F. fulva. These results demonstrate that Cf-4 and Cf-5 share conserved immune initiation mechanisms, but diverge in downstream signaling dynamics, and that the intensity and timing of the HR-related cell death do not affect the robustness of the resistance.

Coumarins are phenylpropanoid-derived specialized metabolites that contribute to plant defence, shape plant-microbe interactions in the rhizosphere, and promote iron acquisition. In Arabidopsis thaliana, a model plant for iron-responsive coumarin metabolism, the enzymatic origin of the catecholic coumarin esculetin has long remained unresolved. Here we identify the first O-demethylation reaction in Arabidopsis specialized metabolism and show that 2-oxoglutarate- and Fe(II)-dependent dioxygenases catalyze scopoletin 6-O-demethylation to form esculetin. We designate these enzymes scopoletin 6-O-demethylases (S6ODs) and validate their activity through biochemical characterization, together with metabolomic profiling and independent loss-of-function mutant lines providing genetic evidence in planta. Disruption of S6OD activity remodels coumarin profiles and alters plant performance under limited iron availability, indicating that esculetin biosynthesis contributes to plant responses under these conditions. Our findings resolve the long-sought missing step in esculetin biosynthesis. It establishes O-demethylation as a previously unrecognized reaction in Arabidopsis specialized metabolism and suggest that 2OGD-mediated O-demethylation is recurrently recruited during evolution of plant metabolism, with implications for metabolic engineering and improvement of iron acquisition traits in crops.

Elongating rice internodes present a developmental gradient from dividing meristem to mature cells, providing an elegant pseudo-time course for study of plant vegetative development. We tested the hypothesis that DNA methylation regulates gene expression during rice internode development by integrating RNA-seq and bisulfite DNA sequencing across eight internode segments. Previously described topologically associated chromatin domain borders aligned with transcription start sites of constitutive expressed genes. CpG and CHG differential methylation was enriched in young segments, consistent with maintenance methylation; whereas CHH methylation showed similar differential abundance in young and old segments. CHH and CHG methylation in upstream regions, CpG methylation within gene bodies, and any methylation in 5' and 3' untranslated regions were permissive of moderate to high gene expression. Very low expression was associated with CpG methylation upstream, CHG and CHH methylation within gene bodies, and CpG and CHG methylation downstream. A nonrandom subset of genes, including cell wall-related glycoside hydrolases, lignin and tricin biosynthesis enzymes, and WD40 proteins, showed methylation-expression correlations, with expression changes enriched in triple-marked elements. These results suggest that internode phenotypes of DNA methylation machinery mutants relate to alteration of specific target genes, opening approaches for grass culm improvement for lodging resistance and biomass production.

Horizontal gene transfer (HGT) allows the movement of DNA across broad evolutionary distances without sexual reproduction. In grasses, HGT is widespread and although a few horizontally transferred genes (HTG) are adaptive, most are purged over time. Within the adaptive HTG, biosynthetic genes encoding enzymes that act together in the same pathway and physically co-localise in clusters have been reported multiple times. The aims of this study are to test whether HGT is bidirectional in a pair of grass species, maize and Zuloagaea bulbosa, and if HTGs are found more than expected by chance in biosynthetic genes organised in clusters. To achieve this, we firstly generated a phased reference genome for Z. bulbosa. Then we identified 56 candidate horizontally transferred genes, of which 45% were from Andropogoneae, including two likely to be of maize origin. Since transfers from Z. bulbosa to maize were previously described, our results show that HGT is bidirectional, although the balance might not be even. After predicting all biosynthetic gene clusters in the Z. bulbosa genome, we found that HTGs are enriched in biosynthetic genes organised in clusters. This correlation between HGT and gene clustering is likely to be a consequence of selection due to the immediate adaptive benefit a whole pathway can provide. Two of the HTGs from Andropogoneae belong to the benzoxazinoid BGC, which previously underwent an ancestral transfer from Panicoideae into Pooideae. The dynamism of biosynthetic gene clusters, including recurrent horizontal gene transfers, contributes to the extraordinary metabolic diversity present in plants. Significance statementHorizontal gene transfer (HGT) is a significant driver of evolution that is widespread in grasses. In this study, we show for the first time reciprocal transfer of DNA between maize and another Mexican grass, Zuloagaea bulbosa. This result represents a proof of concept of bidirectional HGT which allows for limited, recurrent gene flow among distant species. Furthermore, we show that horizontally transferred genes are enriched for biosynthetic genes organised in biosynthetic gene clusters, regions of the genome that encode for multiple enzymes that act in the same biosynthetic pathway. We hypothesise that the transfer of a complete multi-genic pathway, ready to be used and potentially offering an evolutionary advantage, might promote a predominant retention of gene clusters in comparison with background HGT.

Clubroot disease, caused by the obligate biotrophic pathogen Plasmodiophora brassicae, is a major threat to canola (Brassica napus) production worldwide. Clubroot-resistant (CR) cultivars remain the most effective disease-management strategy, but the genetic basis of resistance in commercial canola remains poorly understood because many resistance sources are proprietary and associated genotypic information is rarely accessible. Although nucleotide-binding leucine-rich repeat (NLR) immune receptors account for most cloned CR genes, no pan-NLRome has incorporated CR lines used in commercial canola breeding. Here, we combined whole-genome sequencing and resistance gene enrichment sequencing (RenSeq) to assemble and annotate the NLR repertoires of five homozygous CR inbred lines (IH1-IH5) used for commercial breeding and displaying contrasting resistance profiles against predominant Canadian P. brassicae pathotypes. We integrated these NLRomes with the susceptible cultivar Westar to construct a comparative pan-NLRome for canola. Across the five CR lines, total NLR content was highly conserved, ranging from 504 to 517 genes, with TIR-NLRs representing the predominant class. C-JID-containing TIR-NLRs accounted for more than 30% of each NLR repertoire, and integrated-domain analysis identified conserved and genotype-specific NLR-IDs, including previously unreported domains in IH4. Pan-NLRome analysis resolved 366 NLR orthogroups (OGs), 60.7% of which were core, and identified resistant-line-enriched OGs absent from Westar as candidate CR-associated loci. Unexpectedly, a homolog of the functionally characterized CR gene, CRa, was detected in five CR lines. Moreover, a homolog of another CR gene, Crr1a, was detected in both resistant and susceptible lines, indicating that the presence/absence of a gene alone does not predict resistance. Instead, structural variation affecting LRR and C-JID regions suggests that allele-level diversity within conserved NLR loci contributes to CR-associated variation, with implications for allele-specific marker development and durable CR deployment.

Legumes include some of the worlds most significant crop species, such as cowpea (Vigna unguiculata), a subsistence crop widely grown in sub-Saharan Africa. Despite their importance, legume crop improvement is hindered by a lack of high-resolution expression data, particularly for reproductive tissues and cell types. Here, we report on VigExp, a tool for visualising cowpea gene expression datasets. We demonstrate its utility across a range of vegetative and reproductive cell types of varieties IT97K-499-35 and IT86D-1010, which exhibit 93.75% protein sequence conservation and are amenable to stable transformation. This includes previously published transcriptomes of vegetative, floral and seed tissues, combined with developmentally staged male and female reproductive tissues. Also integrated are novel transcriptomes of laser-captured cell types covering reproductive development from meiosis to early embryo formation post-fertilisation. Spatial expression patterns and transcript levels can be visualised through an electronic fluorescent pictograph (eFP) browser. Validated by RT-qPCR, in situ hybridisation, transgenic, and CRISPR gene editing analyses, the predictive accuracy of VigExp matches prior cowpea functional study observations. Critical genes for nodule development and regulation were also identified and their expression patterns established in cowpea. Novel reference genes, constitutively expressed gene promoters for visualization makers/gene-editing, and tissue and cell specific gene promoters for targeting these regions, are identified. The A-type cyclin, VuTAM2, was also identified, with a critical role in male meiosis established. Collectively, VigExp represents an adaptable and updatable resource to support crop improvement in cowpea and other legumes, which are often highly syntenic with respect to genome composition.

Aquaporins at the plant plasmalemma are divided into two highly conserved subclasses, PLASMA MEMBRANE INTINSIC PROTEINs 1 (PIP1) and PIP2. Arabidopsis thaliana encodes five PIP1 and eight PIP2 isoforms. Individual loss-of-function mutants had been employed for functional analyses. Here, we observe that the pip2;1 pip2;2 pip2;4 pip2;6 pip2;7 quintuple mutant defective of major PIP2 isoforms concomitantly leads to a strongly reduced PIP1 protein level. Lower order mutants pip2;1 pip2;2 and pip2;1 pip2;2 pip2;7 still harbor only 60% and 20% residual PIP1, respectively. This repression is established post-translationally, since neither PIP1s steady-state transcripts nor polysome-associated PIP1 mRNAs are suppressed by pip2;1 pip2;2 pip2;7. Thus, the two major pathways operating in eukaryotes for removal of aberrant proteins, ubiquitin proteasome system (UPS)-dependent ER-associated degradation (ERAD) and autophagy/vacuole-linked degradation, were assessed. Introgression of atg7 blocking autophagy-mediated degradation does not affect the PIP1 protein level of pip2;1 pip2;2 pip2;7. In contrast, introgression of ERAD loss-of-function mutations hrd1A hrd1B and dln1 into pip2;1 pip2;2 pip2;7 partially stabilizes its PIP1 protein level. PIP1 accumulates intracellularly upon pharmacological inhibition of proteasomal degradation by MG132. Nevertheless, the lack of a full PIP1 recovery by these means suggests the flexible operation of parallel ERAD components or unknown pathways. In conclusion, the essential dependence of PIP1 expression on PIP2 isoforms intrinsically interconnects the two PIP subclades at the protein level and will thereby affect their mutual functions. Significance statementPlasma membrane intrinsic proteins constituting the most homogenous plant aquaporin family are nonetheless split into two highly conserved subfamilies, PIP1 and PIP2. The loss of major Arabidopsis PIP2 isoforms does not lead to compensation by PIP1 members, but rather to PIP1s concomitant, post-translational repression. This dependence of PIP1 isoforms inevitably ties the two PIP subfamilies and their function.

Poplar (Populus spp.) is a model system for tree biology. Specifically, P. tremula x P. alba INRA 717-1B4 (hereafter "poplar 717") has become an important platform for functional genomics and synthetic biology due to its rapid growth and ease of transgenesis. Here, we present a single-cell RNA-seq atlas of the poplar 717 shoot, including apical meristem, primary and secondary stems, and three stages of leaf development. Analysis of ca. 159,000 cells resolved 40 transcriptionally distinct clusters representing 7 major cell types, providing a high-resolution view of shoot development and tissue organization. We focused on the epidermis which constituted >15% of cells in the shoot atlas for in-depth characterization of epidermal heterogeneity. By integrating known marker genes with transcriptomic signatures consistent with established poplar leaf phytochemistry, we annotated epidermal cell subclusters corresponding to developmental stages, spatial location, and specialized cell types, including a distinct population of non-glandular trichomes. Coupling the single-cell RNA-seq atlas with bulk transcriptome data from glabrous mutants enabled the identification of novel trichome markers. Experimental validation of a representative trichome-specific promoter established a tool with potential to support cell type-targeted-metabolic engineering. We provide the poplar 717 atlas to the community through the BioPoplar Atlas Viewer (http://bio-poplar-atlas.com), providing a platform to explore the poplar transcriptome at single-cell resolution and a foundation for data-driven cell type-aware genetic engineering strategies in poplar.

Aegilops species are the closest wild relatives of wheat and an important reservoir of genetic diversity for its improvement. Despite their potential, many Aegilops genomes remain poorly characterized. Here we present high-quality assemblies of 18 diploid, tetraploid, and hexaploid Aegilops genomes, which, along with the previously published genomes, complete the production of reference assemblies for all 25 genomes in this genus. Assembly sizes ranged from 5.24 Gb in diploids to 12.65 Gb in hexaploids, with scaffold N50 values up to 749.2 Mb. Gene annotation identified 53,035-156,779 protein-coding genes, of which 21,865-60,490 were classified as high-confidence. Orthogroup-based pangenome analysis across the 25 Aegilops genomes identified 80,521 orthogroups, including 15,809 core, 61,735 dispensable, and 2,977 species-specific orthogroups, highlighting substantial gene content variation among genomes. Phylogenetic analysis of 63 Triticum and Aegilops genomes/subgenomes based on near single-copy orthologs defines the phylogenetic relationships within the Triticum/Aegilops complex and confirms diploid progenitors of polyploid lineages. Ae. mutica (T) and Ae. speltoides (S) belong to the B lineage while the remaining Sitopsis grouped within the D lineage. Structural variation analyses using diploid progenitors as references revealed extensive large-scale rearrangements following polyploidization, emphasizing the dynamics of their evolution. Transposable element (TE) annotation further highlighted subgenome-specific TE expansions and contractions, providing insights into the mechanisms shaping genome structure after polyploidization. Collectively, these genomic resources provide a comprehensive framework for exploring Aegilops diversity, understanding polyploid evolution, and accelerating wheat improvement.

BackgroundLeaves maintain hydraulic homeostasis during photosynthesis through the coordinated action of stomata, which regulate gas exchange and transpiration, and veins, which supply water to the leaf lamina. While functional links between stomatal and vascular traits are known in dicots, their potential genetic coordination in C4 crops remains poorly understood. We investigated the genetic architecture of these traits in maize using a Multi-parent Advanced Generation Inter-Cross (MAGIC) population and a low-cost, high-throughput phenotyping platform integrating leaf clearing, digital microscopy, artificial intelligence, and image analysis ResultsWe phenotyped 285 recombinant inbred lines and the MAGIC founder lines, generating 8,072 images from 2,026 leaf samples taken from seedlings grown in controlled conditions. A YOLOv8-based model automatically detected stomata, while a custom and efficient image-processing pipeline quantified vein traits and stomatal spatial distribution patterns along cell bundles. This enabled simultaneous characterization of stomatal density, size, and distribution together with vein density, thickness, and bundle-associated spatial patterning. Substantial phenotypic variation was observed among genotypes, with strong correlations between abaxial and adaxial traits but no significant correlations between stomatal and vein traits. QTL mapping identified 37 genomic regions associated with stomatal and vein traits, including loci containing known developmental regulators such as stomatal density and distribution1 and stomagen1, as well as novel loci controlling stomatal spatial patterns, divergence between leaf surfaces and veins traits. ConclusionsThese results support independent genetic control of stomata and veins and decoupled contribution to water-use efficiency, providing a novel genetic framework to independently optimize leaf hydraulic capacity and gas exchange in target environments.

Plant regeneration is a major determinant of transformation and genome-editing efficiency, yet the endogenous regulatory networks controlling regenerative competence in horticultural crops remain incompletely understood. The miR319-TCP module regulates multiple developmental processes in plants, but its function in lettuce regeneration has not been defined. Here, we performed a genome-wide analysis of the TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) gene family in lettuce (Lactuca sativa). Thirty-three LsTCP genes were identified and classified into Class I/PCF, Class II/CIN, and Class II/CYC/TB1 groups. Five CIN-class genes, LsTCP2, LsTCP3, LsTCP4, LsTCP10, and LsTCP24, were predicted as high-confidence miR319 targets and supported by degradome-based cleavage evidence. MIR319-overexpression (OX319) explants showed enhanced de novo shoot regeneration, with 94.5% regeneration efficiency and 1.92 shoots per explant, whereas STTM-miR319 suppression (S319) explants showed reduced regeneration, with 28.5% regeneration efficiency and 0.36 shoots per explant. These phenotypes were associated with altered expression of several miR319-targeted CIN-TCP genes, particularly LsTCP4, LsTCP10, and LsTCP24. Disruption of LsTCP4 increased regeneration efficiency to 91.4% and shoot production to 2.05 shoots per explant, resembling the regeneration-enhancing effect of miR319 overexpression. In contrast, disruption of the non-target CIN gene LsTCP17 did not significantly affect regeneration under the tested conditions. Together, these results identify LsTCP4 as a key miR319-responsive negative regulator of de novo shoot regeneration and highlight miR319-mediated repression of LsTCP4 as a potential endogenous strategy for improving lettuce regeneration.

Copper (Cu) is an essential micronutrient whose bioavailability is strongly affected by soil physicochemical properties. During evolution, plants have developed mechanisms to flexibly adjust their metabolism to Cu status. Superoxide dismutases (SODs), including Cu/ZnSOD1 (CSD1) and FeSOD1 (FSD1), are key antioxidant enzymes regulated in Cu dependent manner in Arabidopsis thaliana. Examination of CSD1 cellular distribution and activity revealed that CSD1 is a nuclear and cytosolic SOD whose abundance and activity respond to Cu availability inversely to FSD1. Combined microscopic and biochemical analyses of Cu-dependent dynamics revealed that, unlike FSD1, CSD1 localization in guard cells (GCs) remains independent of Cu availability. CSD1 escapes miR398-mediated regulation in GCs through a cell type-specific splice variant (CSD1.2) that carries an altered miR398-binding site. In silico analyses indicate that this mechanism is also present in crop species. Functionally, the csd1 mutant showed reduced sensitivity to abscisic acid (ABA)-induced stomatal closure, a phenotype rescued by reintroducing CSD1. Biochemical and reactive oxygen species (ROS) level analyses indicate that CSD1.2 most likely acts independently of its canonical enzymatic activity in GCs and functions upstream of the ROS burst in the ABA signaling pathway. Together, we present a novel, cell-type-specific mechanism that safeguards ABA-driven stomatal closure under fluctuating Cu supply.

Flax (Linum usitatissimum L.) has been domesticated for dual end uses as linseed and fiber flax, yet the genomic basis of morphotype divergence remains unclear. Here, we constructed a morphotype-resolved pangenome by integrating three newly generated near telomere-to-telomere genome assemblies with 14 previously published ones. Despite substantial variation in assembly size, driven primarily by DNA transposons, gene content was highly conserved, with little evidence for significant morphotype-specific gene presence-absence variation. Population genomic analyses of 407 accessions revealed that fiber flax had reduced nucleotide diversity, extended linkage disequilibrium, and a more compact population structure relative to linseed, consistent with stronger selection and a narrower genetic base. Genome-wide differentiation was heterogeneous and concentrated in discrete regions. Integration of FST, nucleotide diversity ratios, Tajimas D, and genome-wide association signals identified morphotype-enriched genomic blocks distributed across the genome. Many candidate regions are primarily supported by directional shifts in nucleotide diversity rather than extreme differentiation, indicating selection on standing genetic variation. Genome-wide association analyses identified 1,712 unique quantitative trait nucleotides (QTNs), with predominantly small effect sizes and strong enrichment in gene-proximal regions, consistent with a polygenic architecture. Overall, fiber flax traits tend to be controlled by fewer loci with moderate-to-large effects, whereas linseed traits exhibit a more diffuse genetic architecture. Patterns of Tajimas D further support non-classical selection dynamics, with predominantly positive values in linseed and localized negative values in fiber flax, consistent with selection on standing genetic variation. Together, our results suggest that flax morphotype divergence is driven primarily by selection on pre-existing allelic variation within a conserved gene repertoire. This study provides a comprehensive framework linking genome structure, population genomics, and trait architecture, and highlights the importance of standing genetic variation as a key resource for flax breeding and improvement.

In Asteraceae species, the capitulum is a compact inflorescence, featuring a characteristic reproductive structure. Despite the identification of a few key regulatory factors, the transcriptome-level information on the developing capitulum remains limited. Here, we applied single-cell and spatial transcriptome sequencing to investigate the developing Gerbera hybridas capitulum during floret differentiation. We obtained a transcriptomics atlas encompassing different stages of the Gerbera capitulum and analyzed the cellular and spatial dynamics of gene expression. Using marker gene expression and GO enrichment of cluster-specific DEGs, we annotated putative cell types and described changes in gene expression across sampled stages, potentially associated with ongoing developmental processes. We detected activity of previously undescribed MADS-box genes and defined their spatial expression patterns. Notably, the MADS-box gene GAGL12 was found to be enriched in the putative capitulum phloem cells. The GAGL12 protein was shown in yeast two-hybrid assays to interact with several other MADS-domain proteins with hypothesized functions in vasculature development, and further detailed in silico analyses supported a candidate role in the development of capitulum vasculature. Altogether, we provide integrative and dynamic transcriptomic insight into capitulum and floret development and lay a basis for future functional studies of the control and development of this intriguing reproductive structure.